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Symmetry classification of two-body central potentials in a two-particle Schrodinger 
equation in terms of contact transformations of the equation has been investigated. 
Explicit calculation has shown that they are of the same four different classes as 
for the point transformations. Thus in this problem contact transformations are not 
ON ■ essentially different from point transformations. We have also obtained the detailed 

"pH | algebraic structures of the corresponding Lie algebras and the functional bases of 

O |i invariants for the transformation groups in all the four classes. 

■ 1 Introduction 

a 

> 



x 

5-H 



The position of contact transformations ||, ||] lies in between the point transforma- 
tions and the Lie-Backlund transformations ||, ||. However, in studying the dynamics of a 
system the group of contact transformations has a very important position. Because of the 
continuity conditions in quantum mechanics of the wavefunction and its space derivatives, 
groups of contact transformations play such important roles in the dynamics of physical 
systems. 

Point transformation groups for any set of differential equations involve transformations 
among the independent space-time variables and the dependent variables. The generators 
of the transformation group thus involve only these variables. The groups of contact 
transformations, on the other hand, involve these variables as well as the gradients of 
the dependent variables. The contact relations connecting the gradients with the original 
variables are thus included in the set of differential equations. If the generators do not 
involve the gradients in an essential manner, then the contact transformation group is not 
essentially different from the point transformation group Pj. 
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In a previous work || we have studied symmetry classification of two-body central po- 
tentials for the point transformation groups of time-dependent Schrodinger equation for 
a two-body system. Here, we have done a similar study of the groups of contact trans- 
formations of the same system. What we have found here, and what could not be known 
without this detailed analysis, is that the groups of contact transformations for this system 
are not essentially different from the corresponding groups of point transformations. Thus 
we have four classes of two-body central potentials: 

1. a constant potential with a 31-parameter Lie group, 

2. a harmonic oscillator potential with a 20-parameter Lie group, 

3. an inverse square potential with a 16-parameter Lie group, and 

4. all other potentials with a 14-parameter Lie group. 

Using Lie-Jacobi's method Q |[ we have calculated the functional bases of invariants 
for all these cases. Any invariant of these groups can be functionally expressed in terms 
of these base invariants. For the constant potential case there are 5 base invariants. The 
corresponding number for all the other cases is 4. 



2 Contact transformations and Schrodinger equation 

The method H] for obtaining the group of contact transformation is a generalization of 
the method of extended group ||, ||] for obtaining that of point transformations. We give 
below the essential points as it appears in our system of 2 particles of masses m\ and ni2 
at positions f\ and r*2- 

We use the relative space coordinate r = r\ — r-i , the centre-of-mass space coordinate 
R = (mifi + m2r<2)l(yn\ + 777,2), the reduced mass m = mxm^/ [m\ + 771-2) and the total 
mass M = m\ + 777,2. It is to be noted that < m/M < 1/4. The limit m/M = 1/4 occurs 
when the two particles have equal masses. Positronium atom and a homonuclear diatomic 
molecule are two important physical systems which have this limiting value of m/M. The 
other limiting value m/M = occurs when one of the masses and hence M is 00. This is 
same as ignoring the motion of the centre-of-mass. This does not describe a true two-body 
system and cannot thus throw any light on the classification of inter-particle potentials. 

In terms of these variables, the 2-particle Schrodinger equation for the wavefunction 
\E r ° becomes 

a A 

with the contact conditions 

A ca = ^° qCCt - ^> ca = 0, (2) 

A rA = ^J rA - * rA = 0. (3) 

We have written the physicist's version of the equation, keeping the Planck's constant, 
the relevant masses and the imaginary number 1. We use a compact notation q a , a = 0, 
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ca, rA, where q° = t, q ca = R a , q rX = r\. Here and later on a, A etc. will mean the 
cartesian components, subscripted variables (other than the cartesian components) will 
mean derivatives with respect to those subscripts, the letters c and r will mean center-of- 
mass and relative coordinates, and t will mean time. 

Equations (||), @ are actually the definitions of the gradient variables. The inter- 
particle potential v(r) has been taken to be of central nature. In terms of the gradient 
variables the Schrodinger equation is now of the first order. 

The generators of the Lie group transformations in the space of q a , *$> a are of the form 



* = E 



(4) 



with a = 0, ca, rA. The arguments of £ a and \ a contain the collection of all q a and *$> a . In 
the Racah nomenclature the £s and x s are called the velocity vectors of the generator X. If 
X° explicitly contains the gradients ^> ca , ^> rX , then the contact transformation is essentially 
different from a point transformation [g]. 

As in the case of point transformations, the first extention of X is written as 

q b 



Here, 



x a;b = «-E ^4 - E KK^' + E ■ ( 6 ) 

b> a',V a' 

The effect of X^ 1 ' on the different As of equations (Jl|), (Q), (Q) are 

A° - ix«* + ^ E X™ + ^ E ^ - «Mx° ~ ^ E ? X r& = °* (7) 

a A A 

I^'A™ = X °'' ca - X ca = 0, (8) 
XW A rX = x °'' rX - X rX = 0. (9) 



3 Defining equations of group generators 

The group of contact transformations for the equation (|l|) will be uniquely known when the 
velocity vectors £ a and \ a win De obtained. These velocity vectors satisfy an overcomplete 
set of differential equations known as the defining equations. These defining equations 
are obtained by separately equating to zero the coefficients of the different monomials in 
^f a b appearing in equations ([?]), (|§|), (||). 
In our case we get the defining relations: 

= = 0, (10) 

+ ^ ca ^o = ej-A + * rA ^ = o, (ii) 



54 



P. Rudra 



V 

A\I> ca 



E* c ^-E* rM ^ = 



/3 M 



x^-E* c %--E* r %- 



x ^-E^-E^ 



x ca = o, 



x^-E^-E^, 



r/i 



r/i 



- x rA = o, 



U Q7 

ft 

2M 
ft 



/3 H 



2m 



X *o-E^fo-E^Co-^ o 

/3 

_ft_ 

2m 

ft 
2M 



XvJirA 



rA^ca 



+ $L - + * ca a = o, 



+ C^-x^ + * rA ^o = o, 



rA 



= o, 
= o, 



x?o-E* c/ ^ c o-E*% r o 



ft 



+ ■u(r)* 
ft 



(13) 
(14) 
(15) 
(16) 

(17) 



(18) 



(19) 



(20) 

(21) 

(22) 
(23) 



x^-E^S-E*^ 



^0 



a A 



-u(r) 



x° + *°& ' V{r) ' •"° 



A 

. ^°E^e 

rv(r) 



(24) 
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We have obtained the general solution of these defining equations and have indicated the 
method in the Appendix 1. The velocity vectors are expressed in terms of the auxilliary 
functions 



F°(t, R, r) = Mt) -jE fo a (t)'Ra - ^ E tfW'rx 

"a A 



a A 

F ca (t, R, r) = fS"(t) - \b{t)'R a + £ e Q ^ 7 /r 7 ^ + m £ / CQ;rA r A , 

/3 7 A 

i^ A M,r) = / rA (t) - i&(*)v A +x) e wr% - M Er%. 



and are of the form 

e(t) = b(t), 

e a (t,R,f) = -F ca (t.R,r), 
e X (t,R,r) = -F rX (t,R,f), 



X °(t, R, r, *°) = f(t, R, f) + *°F°(t, R, f), 



^ = dx°(t,R,r,V ) +m y f c Q/ «;r, 



<9r A 



where F°(t,R,r) and f°(t,R,f) satisfy 



dF°(t,R, 

1 dT 



h ^ d 2 F°(t,R,r) J^ s ^d 2 F°(t,R, 



2M 



2m (3rA) J 

A 



-t,( r ) & (t)' + ^ e r ^ prX ^ £ f ) = °> 

A 

. df°(t,R,r) h ^ d 2 f°(t,R,r) H ^ d 2 f°(t,R, 
% 8t 2M ^ (8R n ) 2 2m ^ (dr x ) 2 

A 



-v(r)f(t,R,r)=0. 



(25) 

(26) 
(27) 

(28) 
(29) 
(30) 
(31) 

(32) 



(33) 



(34) 



(35) 
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]\[q m ci / 
IN cLlllti / 

Symbol of the generator 


Form of the generator 


Scaling : 


£hT<0 / ; £>\lrca / y £>\lrrX 
a ' ' X 


time translation : 
X 1 


d 


centre of mass coordinate 
space translations : 


. d 

l dR a 


relative coordinate 
space translations : 

yrX 

JVrp 


. d 

dr x 


centre of mass coordinate 
Galilean transformations : 

yea 
X G 


"? + T 


[« + S °9J»] 




relative coordinate 
Galilean transformations : 

yrX 
X G 


tx T +- 


A S ^ d* rA 




centre of mass coordinate 
space rotations : 

X R 


^ ^ ^ C-a(3~{ 
Pi 




relative coordinate 
space rotations : 


-i^ex^v 


_ M <9rv 9*™ 




cross-rotations : 


■ \ ® jca ® 

l yM[ rx dR a d^ r \ 


. [M\ d x d 1 
+ Z V to [ Q 9r A J 


X! 


2tX l - Y R a X r T a - Y r * X T + IvotXs + i*°-^o 

a X 




t 2 X l ~ 

— 1 


tYn<* x T -tY r * x T X 

a X 
a X 

T o •+ ® M \ - D <9 to 

a 


V 9 


relative vibrations : 

yrX 

A v,(±) 









Table 1. Generators that describe the different Lie Algebras of the four classes of inter-particle 
potential. 



Here, a prime on the function of a single variable denotes derivative with respect to that 
variable, e a/ g 7 and e\^ u are the permutation symbols and ff a , /{ A , /g Q ' rA are constants. 
Since x (^-^;^^°) does not contain and ^ rX , the group of contact transformations 
for this system is not essentially different from that of point transformations. However, 
the generators now contain derivatives with respect to the gradient variables, because the 
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group of contact transformations is the first extension of the group of point transforma- 
tions. 

Equation (35) is nothing but the original Schrodinger equation and f°(t,R,r) is a 
solution of equation ([!]). This is the symmetry corresponding to linear superposition 
principle of the Schrodinger equation and forms an infinite dimensional invariant subgroup 
of the total group. The factor group modulo this subgroup is the physical group of interest 
and will be referred to as the group of contact transformations. 

Substituting F°(t,R,f) from equation ( p5|) in equation ( |34|) and equating coefficients 
of different monomials in the space coordinates to zero, we get 



b(ty 



o, 



(36) 



^f \t)" + ^-fS\t) = o, 



v(r) + —rv(r)' 



b(ty = o, 



(37) 
(38) 

(39) 



and either 

v(r)' = 



or 



Jo — u - 



(40) 
(41) 



4 Symmetry classification of interparticle potentials 

^From the solution of equations (|36|)-([4"l|) we get four different classes of interparticle 
potentials. This complete symmetry analysis of the 2-particle Schrodinger equation as 
far as the dynamics of the system is concerned shows that contact transformations do not 
enforce any more restriction than that already required on the basis of point transformation 
symmetry of the system. Their group algebras are the first extensions of the algebras of 
point transformation symmetries of the corresponding potentials || . 

These group algebras are described in terms of the generators given in Table 1. 

The letters T, G, R and V in the symbols of the generators denote space translational, 
galilean, space rotational and vibrational modes described by these generators. We have 
kept the imaginary number i in the forms of the generators so that these can be identified 
with the usual quantum mechanical operators for energy, linear and angular momenta. 

In actual calculations with the group algebras, the structure constants (commutation 
relations in physicists' parlance) are of the greatest help. They are given in Appendix 2. 

Any Lie Algebra L is a semi-direct product of its Radical R and a semisimple part L/R. 
The semisimple part L/R is again a direct sum of ideals which, as subalgebras, are simple 
0, [|. These characteristics of the Lie Algebras for the different classes of inter-particle 
potentials are given in Table 2. 

It is to be noted that the generators for the space translational and Galilean as well as 
vibrational transformations always belong to the Radical of the Lie Algebra. 
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Inter-particle 
potential: 


Constant 


Harmonic 
oscillator 


Inverse square 


Arbitrary 


v(r) 


= v 


2 2 
muj r 

= V0+ 2h 


Vl 

= V - 


7^ "0, 








2 2 

muj r 

Vl 

Vo 2 


Generators of 
Lie Algebra, 
L 


-(As, A ,A GTi j, 
yrX vcoi\r\ 
G,T,R' ' 

Xi,X 2 ) 


f v" yea 
|As,A , A Gril , 

'^v,(±)j 


{As,A*,A G a Ti? , 
Xr X , Xi,X 2 } 


{A s , A*, 

yea V r ^\ 


Dimension 


31 


20 


16 


14 


Solvability, 
Nilpotency, 
Simplicity, 
Semisimplicity 


AT„ 

IN one 


AT,, „ 

IN one 


AT„ 

IN one 


IN one 


Centre, Z(L) 


{Xs} 




{Xs} 




Radical, R 


f V yea yrA 1 
i A S, A G,T> A G,TJ 


f V yea yrA 1 

v*-S)-*g,t> a v,(±)j 


{As, A G a T } 


{As, A G Q T } 


Semisimple 

part as 
direct sum 

L/R = E®h 

of simple 
ideals 


A = {-X" ) -Xi) -X2}, 

r f yea yrA 


= {x*}, 

h = {Xf?}, 

h = {xg} 


h = {X 1 , X\, X2}, 

h = {x%>}, 
h = {x r £} 


A = {x*}, 
h = {x R a }, 

h = {X$} 


Cart an 
subalgebra 
H 


{Xs,X R ,X R , 

X c3;r3 ,Ai} 


{Xs,X l ,X^ TR , 

x r R 3 } 


{X s ,X^,X r3 , 
Xi}, 


{As, A*, 

yc3 V r 3\ 



Table 2. Characteristics of the Lie Algebras of the Symmetry groups of the different classes of 
inter-particle potential. 



We note that the Radical, being solvable, has only 1-dimensional irreducible representa- 
tions (irreps) pj. Thus the irreps of L are obtained if the irreps of L/R are known. To this 
end we give in Table 3 the algebraic characteristics of the simple subalgebras appearing 
in the direct sum of L/R for the different classes of inter-particle potential. 

5 Functional bases of invariants 

In order to investigate integrability of a dynamical system we require the functional bases 
of invariants in terms of which all invariants of the dynamical system can be expressed 
functionally. These base invariants of a Lie algebra L generate the Centre of the Universal 
Enveloping algebra of L. Since they commute with all the generators of L, they appear as 
conserved quantities of the system. Lie's method |], [H], |ll], O has been utilized to obtain 
the base invariants for the four symmetry groups of Section ||. If the functional base has 
s invariants Ji, I2, ■ ■ ■ , I s , each a function of the r generators X a of the symmetry group, 
then 

[X a ,I b }=0, a = l,...,r, 6 = 1,..., a. 
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Algebra: L 


Li 


Li 


is 


Li 


Generators: 
{Xi} 


{X\X U X 2 } 






m 


Dimension 


3 


15 


3 


1 


Cartan 










subalgebra: 


{-H 


{X*X$, X^} 


{Xr} 


{x*> 


11 










Rank 


i 


3 


1 


1 


-Daae ui 
simple 
roots : 
A 


a = 1 


ai= (i) ,a2= ( ?)' 

V / 


a = 1 


void 










Root 
system 


±a 


±ai, ±a 2 , ±a 3 , 
±(ai + a 3 ),±(a 2 + a 3 ) , 
± (ai + «2 + 03) 




void 


Isomorphy 
to 


0(3) 


0(6) 


0(3) 


U(l) 


Standard 

sctof 
generators 


E a = X 2 
E- a = X 1 


+ (X cl ^ 3 T *X c2 ^ 3 ) , 

£±a 2 = {XgTiXg) 

- (X cl ' r3 T iX c2 ^) , 
E± a3 = (X cl ' rl - X c2 ' r2 ) 

±i (X cl;r2 + X c2 ' rl ) , 

- (x c3 ; rl ± *x c3 ^ 2 ) , 

E±(a 2 +a 3 ) = ( X R ± *^]? 2 ) 

+ (X c3 - rl ± *X c3 ^ 2 ) , 

^±(ai+a 2 +a 3 ) = (A cl ' rl + A c2 ' r2 ) 




void 



Table 3. Characteristics of the simple subalgebras forming the semisimple parts in Table 2. 



If I is any other invariant so that 

[X a ,I]=0, o=l,...,r, 
then / can be functionally expressed as 

I = I(h,...,I s ). 

If the base has the constant as its only member, then the system is completely chaotic. If 
on the other hand s = r and all the generators give mutually commuting invariants, then 
the system is fully integrable. Actual dynamical systems are almost always in between 
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these two extreme cases. In the four symmetry classes obtained in Section f|, the same 
thing happens. 

In Table 4 we define auxilliary operators in terms of which the functional bases of 
invariants for the four classes of distinct inter-particle potentials are given in Table 5. 



Inter — particle 
potential 


Auxilliary operators 


Constant 
potential 


Y R = X S X C R + - V e a/37 l^lg 7 , 

lVl 

h 

Yr X = X s X r R + — ^2 ex^X^X^, 

Y X S X + [ X G X T A T A G J ) 

y* = x s x> + VQ (x s f + A J2 (x?f + A £ (a-) 2 , 

\\ = X s X t - ii (X s f + -^Y1 X " x o + ^ Y. X T X a, 

a A 


Harmonic 
oscillator 
potential 


Y R a = X s X c R a + — ^ ea^^r^G' 

01 

Yr - X S X R + 2 ^ e ^ X V,(+) X v,{-y 

Yt = X sX> + A £ (AT)2 + ±_ ^ X^X^ 

a A 


Inverse 
square 
potential 


Y R a = X s X c R a + — ^ ea^X^X^, 

01 

Y* = A S A 4 + v (X s f + A (*t) 2 , 

a 

Y 1 = XsX, - 4* (X S ) 2 + A £ AT AT, 

a 

f 2 = x s a 2 + A^ (x -)2 

a 


Arbitrary 
potential 


^ft* = ^sA£? + — Y eapyXtfxg, 

01 

Y t = X 8 X t + A^(x-) 2 

a 



Table 4. Auxilliary operators for the four distinct classes of inter-particle potentials. 



In all these cases Is is the scaling operator and I* is essentially the energy operator, 1^ 
and are the cm. coordinate and the relative coordinate angular momentum operators. 
Except for the case of constant potential, these are the only invariants. The only potential 
that has more invariants is the constant potential and the two extra invariants Jg' and 
of degrees 3 and 4 in the generators make this apparently simple system actually more 
complicated yet more systematic. 
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Inter — particle 
potential 



Constant 
potential 



Harmonic 
oscillator 
potential 



Inverse 
square 
potential 



Arbitrary 
potential 



Base Invariants 



Is = Xs, 

I 1 = Y 2 Y t — i (Yi) 2 , 

i° R = J2 «n 2 + E ( Y n) 2 + E ( YCa:rX ) 2 

a X aX 

jc\r \ A \^cayrca;rXyrrX 

3 / j R R 

~\ E E e a (3 7 e x ^Y ca ' rX Y c ^Y c ^, 



af3j Xfiv 



E( r «") 2 



EW 



+E 

a 



ca;rX\rrX 



E 



Y^Y^Y 



ca;rX 



aX 



1 2 



01 v* v 



E y ; 



n 



0i 



Is = Xs, 
jt = Y \ 



E« 

a 

Eft 



D 2 - 



= Xs, 
I t =Y 2 Y t -](Y 1 f, 



E« 



^ = E W 



Is = Xs, 

jt = y t) 

^ = E« 



ca\2 



EW) 2 



Table 5. Functional base of invariants for the different classes of inter-particle potentials. 

It is to be noted that in all the expressions in this section involving products of operators 
symmetrized forms have to be taken. 

Acknowledgement: I thank the referee for suggesting a number of improvements to 
the manuscript. 
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Appendix 1 

In this Appendix we indicate how the defining equations (|l0|)-(|2"4|) are solved to obtain 
the group generators. 



Equations @, (|TTJ) give £° = £°(t). From equations it follows that 

i ca = £ ca (t, R, r, f' A = C X (t, R, r, *°). (42) 

Equations ([R]), ( |T7[ ) thus give 

X°^X°M,r,*I/ ). (43) 

Using the auxilliary function 

F(q a , ^ a ) = x° ~ ~ Yl ^ rX C X , (44) 



A 



equations (|i"8|) and ( |i"9|) now become 



OF 

&F OF 

x°=^-E^^-E^ A ^- (^) 

« A 

Here F is linear in \& CQ and Vl/ r \ Because of this linearity, equations (ll~8|)-(|2l|) give 

r = e a it,R,r), e x = e rA M^, (50) 

and equations (p4[), (H^), ( |48|) and ( |50| ) imply linearity of x° a and x rX with respect to ^ ca 
and \E' rA . Equations fl20|), fl21|) together with equations (fi4|)-(|4"8|) now give 

9 q (51) 

and equations fl22]), ( ]23] ) become 

mC q L + = 0. (52) 



From equation (|2j) equating separately to zero terms quadratic in, linear in and inde- 
pendent of ty ca and ty rX , we get 



0, (53) 
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K q 



ix° q0 



\ " tea 

2M ^ V^s' 



& tea , A v 



2m 



2M 



+u(r)* c 



These will result in 



Lo ~ X*o H ^-r > r A 4 

y (rj 

A 



£° = b(t), e a = -F ca (t, R, f), C rX = ~F rX (t, R, r) 
X = f (t,R,7) + * F°(t,R,r), 



(54) 
(55) 



(56) 



(57) 



with 



F Z + ^ Q 



-b(t)'5 a p, 



F^ + F^ = -b(t)'5 XlM , 
MF™ + m F^ = 0, 



and 



c« 

r u ,r u i 



2m ^ 



A 



- ~^T7 £ /iL-Rc, - 7T— £ /r A ,r A + ^(0/°- 



(58) 
(59) 
(60) 

(61) 

(62) 

(63) 

(64) 



Equations (63) and (|64[ ) are obtained from the terms that are linear in and independent 
of in equation (pq). Equation (58) gives us the result that F ca is linear in R; and 



equation (|5^) shows that F rX is linear in r. Equation ( |60| ) will now give F°, F ca and F rX 
as in equations (|25|)-(|2"7|) and x ca , \ rX as m equations (0), (^), the difference is that 
/f 7 > f[ v and /o"' r are yet functions of time t. Equations (^|), fl62] ) now ensure that these 
are indeed independent of t. 

From equation (|63|) we also obtain equations (|36])-(pT|). 
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Appendix 2 

In this Appendix we give the structure constants of the generators given in Table 1. 



[x t ,x%]=iXF, 



[X t ,X 2 ] =iX 1 +AX s , 



X\X r G X 



iX r T x , =2i(X t + v X s ) , 



Xca vc/3 



iM S Y 



T ; A R 



: Y, ea^xtp, [xp, x,\ = ixp, [X%\X 2 ] = iXg 



X^,X cf3]rX 



iJ^5 a/3 X r T \ 
V m 



-yea -yep 

G > A R 



[Xg* , Xi] = -iXg* , [xg* , X c/3;r A 




- ni—SapXg, 



rn 



A i? ' A i? 



vrX sr r V- 
A T , A G 



%m R Y 



vrX -v-r/j. 
A T ,A R 



Ay , Ai — !Aj , Ay , A2 



7'A 



A G ; A R 



A G 5 A l 



-iX r G \ 



^A/i^ , 



x rX^ x ca;rn 



A^\ A^l = » ^ e^Xfc", fe\ X-H = i Y ex,** 



X ca;r \ X c ^] = il S a pY ^ X R + <V E e ^f X R 

\ v 7 / 

[X 1 ,X 2 ] = 2iX 2 , [x\X r v x {±) 

Xr , Xy^ = I E e XllV X V,{±) ' 



TuX rX 



v,(±)> 



vrA 

A ^(+)> A y,(-) 



2muj 



Sx/j,Xs. 
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